The split-ring Josephson resonator as an artificial atom 
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Using the resistive-shunted-junction model we show that a split-ring Josephson oscillator or radio- 
frequency SQUID in the hysteretic regime is similar to an atomic system. It has a number of 
stationary states that we characterize. Applying a short magnetic pulse we switch the system from 
one state to another. These states can be detected via the reflection of a small amplitude signal 
forming the base of a new spectroscopy. 



PACS numbers: Josephson devices, 85.25.Cp, Metamaterials 81.05.Xj, Microwave radiation receivers and 
detectors, 07.57.Kp 
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I. INTRODUCTION 



Since the realization of lasers using atoms in the six- 
ties, many researchers have been trying to replicate this 
effect using quantum objects similar to two-level systems. 
The basic elements to realize a laser are [l| (i) energy lev- 
els for the electrons, (ii) a pumping mechanism to pop- 
ulate upper levels and (iii) an optical cavity to confine 
photons enhancing the probability of interacting with 
excited electrons. A Josephson junction between two 
superconductors is a macroscopic quantum system and 
Tilley predicted that interconnected Josephson junc- 
tions could emit a coherent radiation. Rogovin and Scul- 
ley also established a connection between Josephson 
junctions and two-level quantum systems. The super- 
radiance prediction was confirmed by Barbara et al f4| 
who showed that the power emitted by an array grows 
like iV^ where N is the number of active junctions in the 
array. See also the recent experiments by Ottaviani et 
al [5| showing the synchronization of junctions in an ar- 
ray. Since the radiation emitted is in the Terahertz range 
where there are no solid-state sources, these systems are 
being investigated as microwave sources (see for example 
0). Up to now however practical difficulties subsist and 
the powers emitted by such devices remain low. 

From another point of view artificial materials (Meta- 
materials) have been fabricated by embedding metal 
split-rings or rods into dielectric materials. This way neg- 
ative index materials have been fabricated^, One 
can obtain artificial atoms with a high magnetic moment 



and a nonlinear response to electromagnetic waves. Such 
an artificial atom would have many advantages over a real 
atom. For atoms the dipole momentum is very small and 
the interactions between them are small. To realize in- 
teractions large density are necessary. Also in general the 
excited state has an energy much larger than the energy 
of the dipole coupled with electromagnetic wave. Many 
candidates for artificial atoms have been proposed, most 
of them width intrinsic nonlinearities. Among them we 
have a diode , a Kerr material 11 1 or a laser am- 
plifier fl3|. Another example is the Josephson junction 
discussed above (see for a review). The use of these 



devices in metamaterials was advocated by Lazarides|14 
and by the authors [l^. They introduced a split 
ring resonator with a Josephson junction contact. In the 
Josephson community, this device is called an RF SQUID 
for Radio Frequency Superconducting Quantum Interfer- 
ence Device [£] . It is the elementary component of the 
arrays of junctions discussed above. In this work we will 
show that rather than the junction itself, the RF SQUID 
can be considered as an artificial atom. 

We will study the so-called hysteretic regime where the 
system has controlled metastable states and show that 
one can switch from the ground state to one of these 
excited states by applying a suitable fiux pulse. We show 
that a magnetic field can act on such a system in a similar 
way as an electric field acts on dipoles in atoms. We also 
show that these states can be detected by examining the 
reflection coefficient of an electromagnetic wave incident 
on the device. This is the base of spectroscopy. The main 
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where V and / are respectively the voltage and current 
across the barrier, are the macroscopic phases in 

the two superconductors, Ic is the critical current of the 
junction and 4>o — fi/(2e) is the reduced flux quantum. 
The right panel of Fig. [T]sliows the equivalent electric cir- 
cuit of the whole system assuming a Resistively Shunted 
Junction model [3],[1| for the Josephson junction. We 
'-then define 



V" 



FIG. 1: A split ring resonator with an embedded Josephson 
junction (left panel). The right panel shows the equivalent 
circuit using the Resistively Shunted Junction model for the 
Josephson junction. 



result of this study is to show that a split-ring Josephson 
oscillator (RF SQUID) in the hysteretic regime behaves 
as an artificial atom with discrete energy levels. It is the 
only device that leads to such discrete levels as opposed 
to the systems mentioned above. 

The article is organized as such. In section 2 we derive 
the model and analyze it in section 3. In section 4 we 
characterize how to switch from one state to another and 
how the state can be detected. In the last section we 
study the scattering of an electromagnetic wave by a split 
ring Josephson resonator 



II. THE MODEL 

The device we consider is shown in the left panel of 
Fig. [T] It is a split ring resonator in which is embedded 
a Josephson junction. Practically it can be made using a 
ring like strip of superconducting material where a small 
region was oxidized to make the junction. The right panel 
of Fig. [1] shows the electric representation of the device, 
an inductance L for the strip and the Resistive Shunted 
Junction (RSJ) model for the Josephson junction. The 
latter represents the junction as a resistor i?, a capacity C 
and the nonlinear element in parallel. This last element 
is the sine coupling Ic sin ^/(j)Q where $ is the magnetic 
fiux and (f>Q is the reduced fiux quantum (see below). In 
standard electronics the conjugate variables are voltages 
and currents while in superconducting electronics they 
are the fiuxes and charges where a flux is defined as the 
time integral of a voltage. This is why we present the 
derivation in detail here. The device is assumed to be 
operating at low temperature so that losses are minimal 
and Josephson relations hold. 

The Josephson equations describing the coupling of 
two superconductors across a thin oxide layer are 



V = — ($' - $") 
at 



I — In sin 



(1) 



T/'(r)dr, $" 



y"(T)dT, 



where V , V are the voltages on each side of the Joseph- 
son junction (see Fig. [1]). Kirchoff's law gives 



V -V' =-LIt +Ee, 



(2) 



where the subscript indicates time derivative and is 
the electromotive force due to an electromagnetic pulse 
incident on the ring. We neglect the resistance of this 
loop which we assume to be made of superconducting 
material. In terms of fluxes this relation is 



-LI + 



(3) 



where $e is the incident flux. Kirchoff law for node V 
gives 



Ic sin( 



Do R 
which in terms of the fluxes becomes 

-) + — 



$" V' - 

) + ^^ + CiV;-V^\), (4) 



I — Ic sin(- 



R 



+ c(<i>;, -$"«)• (5) 



We now introduce the phase difference $ = $' — and 
combine equations ([3]) and (O to obtain our final equa- 
tion 



— +IcSm{ — ) 
R 00 



(6) 



The quantity in brackets is the current / circulating in 
the loop. To measure the importance of the sin term in 
this equation we introduce the dimension-less Josephson 
length as a ratio of the flux in the loop versus the flux 
quantum 

00 

Time is normalized by the Thompson frequency, t' = LOTt 
where 

UJT = -^=- (8) 



The fluxes are normalized by 0o as $ = 0(/>Oj = 4'e4'o 
In the normalized time t' , dropping the ' for ease of writ- 
ing we get our final dimensionless equation 



(j)tt + Oi4't + P sin(0) -I- = 0e, 

where the damping parameter a is 

lutL 



(9) 



(10) 
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III. ANALYSIS OF THE MODEL 

The ordinary differential equation (|9]) can be written 
as the 1st order system 

<f>t = ^, (11) 
V't = -aV'-/3sin((/))-0 + (/)e. (12) 

The system has the fixed points (0,0) and (0*,O) where 

- /3sin(</>*) - 0* +(/)e = 0. (13) 

A plot of the above relation indicates that for /3 > 4.34 
there are no additional fixed points. The fixed points 
can be approximated for large /? using an asymptotic 
expansion. The equation ()13p can be written as 

sin0* + i(^*-0,) = O. 

Writing the solution 

1 l2 

(j>* = (1)0 + -(j)i + - (1)2 + ■■■ 

we get 

(f)* ^ nn + ^i-iy\,p, - mr) + . . . , (14) 

where n is an integer. 

In the absence of damping a — and forcing (/>e = 0, 
the system is Hamiltonian with 

F(0,0t) = ^0?+/3(l-cos0) + i02. (15) 

The stable fixed points correspond to the minima of the 
potential 

F(0)=/3(l-cos0) + i02, (^g) 

Fig. [5] shows a plot of the potential V{(t)) for /3 = 1, 9.76 
and 100. For (3 = 1 shown as a continuous curve (red 
online) there is only one fixed point = 0. For /? — 
9.76 shown in dashed line (green online) there are three 
minima corresponding to stable fixed points, <j> — 0, ±0* 
where (j>* « 2tt. For [3 = 100 there are many stable fixed 
points. 

Another point is that the incident flux can be used to 
modify the energy levels of the system, assuming the 
incident flux to be constant we can add a term to the 
potential and obtain the generalized potential 

1/(0) =/?(l-COS0) + (17) 

where 0e is the incident flux, assumed constant. This 
expression is plotted in Fig. |3]for /3 = 15 and 0e = 0, I.Stt 
and 4.57r. The minima are symmetric for (j)^ = and they 
are shifted to the left and the corresponding value of the 




FIG. 2: Potential energy y(0) = /?(l-cos<;/)) + i0^ with three 
different values of /3, /3 = 1, 9.76 and 100. 
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FIG. 3: Potential energy y(0) = /3(1 - cos^) + + h(f> 
for three different values of the static incident flux <j!>e = in 
continuous line (red online), 0e = I.Stt in dashed line (green 
online) and <j)e = 4.57r in short dashed line (blue online). 

potential is decreased. By applying a sufficiently large 
continuous field one can then shift the system from one 
state to the other. 

For this one degree of freedom Hamiltonian, the orbits 
are the contour levels of the Hamiltonian. An important 
orbit is the separatrix connecting the two unstable fixed 
points 0* « TT. It is given by 

i0? + /3(l-cos0) + i02^i/(r,O). (18) 

This value of the Hamiltonian can be approximated for 
(3 » 1 as 

i/(0*,O)«y(l + |)+/3(l + cosp«^ + 2/3+^. 

(19) 

Fig. |4] shows the phase portrait for /3 — 9.76. For 
this value there are only five fixed points. Notice the 
closed orbits around the fixed points, the closed orbits 
surrounding the three stable fixed points. 

Wc have shown that the steady states of the split-ring 
Josephson oscillator are similar to the stationary states 
of atoms. The values V{(j)) near 2n7r are the analog of 
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FIG. 4: Phase portrait (<j!>, (\>t) of the Hamiltonian system (fT5| 
(a = 0) for /3 = 9.76. The contour levels presented are 0.1, 1 
5, 17., 25.0360499136927 (separatrix) , 30.,40. and 50. 



atomic energy levels. In the quantum regime we should 
observe Metastability of these states, i.e. there should 
be quantum tunneling through the barriers at the rni 
positions. This would result in a finite life time of these 
steady states. In the next section, we will select the in- 
cident flux to move the system from one equilibrium 
to another. 



IV. INFLUENCE OF DAMPING 

We consider now that the state of the system can be 
shifted from one fixed point to another via an incident 
flux. For a short lived perturbation, the system then 
relaxes freely to a minimum of energy. The influence of 
the damping is essential , it should be present to allow 
the relaxation but small to preserve the picture of the 
potential. To examine how an incident flux will shift 
the system from one equilibrium position to another it 
is useful to analyze the work equation. To obtain it, we 
multiply © by 0t and integrate over time. We get the 
difference in energy 



E{t^)-E(t{) = [-<^?+/3(l 



-COS0) + -0"|;: = / dt(j)e(j)t- 

(20) 

The first term on the right hand side is the forcing while 
the second one is the damping term. When a square 
pulse is applied to the system, such that 

(j>e{t) ~ a, for tl < t < t2, elsewhere 

the first integral is a[(/)(t2)— If the system is started 
at (0, 0) in phase space so that the E{ti) = and (j){ti) = 



E{t2) = a(l>{t2) - a / di</)2, 
Jti 

so that 4'{t2) determines how much energy is fed into 
the system. When the pulse is long t2 >> ti (f>{t) will 
relax and oscillate so that there are values of t2 such that 
0(^2) 0. In that case no energy gets fed into the system. 
A sure way to avoid this is to take a narrow pulse. 

The natural frequency of the oscillator around the 
(0, 0) fixed point is 



Wo 



(21) 



which for j3 = 9.76 gives ujq « 3.28 and a period 
To = 27r/ajo ~ 1.91. To simplify matters we now con- 
sider a pulse of duration much smaller than Tq. This is 
experimentally feasible and can be modeled using a Dirac 
delta function (j)eit) = a6(t), where a is a parameter. Let 
us analyze briefly the solution. The equation ([9]) becomes 



(f>tt + (a + 6)(l)t + 13 sm{(f>) + (f> = aS{t). 



(22) 



Integrating the equation on a small interval of size e 
around 0, we get 

[0t]le + aMle+ / dt(/?sin(/) + 0) = a. (23) 



We now take the limit e — ^ 0. We will assume continuity 
of the phase so that [0]!^ 0. The third term being 
the integral of a continuous function tends to when 
the bounds tend to 0. Assuming (/)t(0_) = we get 
^t(0+) = a so that such a short incident pulse will just 
give momentum to the system. 

We will now explore systematically the plane (a, a) 
characteristic of the incident pulse. The equation © has 
been solved numerically using a Runge-Kutta algorithm 
with step correction of order 4 and 5. 

The plot in the (a, a) parameter plane shown in Fig. [S] 
shows the final states, O the central focus {+), R the right 
focus ( X ) and L the left focus * reached by the system. 
Notice how these are organized in "tongues" following 
the sequence OROLOROL... as one sweeps the 
plane counterclock-wise starting from the horizontal axis. 
This simple geometrical picture can be understood by 
examining Fig. 21 In the case of small damping, the 
iQDaratrices around the fixed points are not affected very 
mui^lt^^heir perimeter is proportional to the probability 
of reaching one fixed point or another. The system is 
moving clock-wise along the orbits. Assume the system 
reaches the central point O for a given set of parameters. 
If the damping is increased, the orbit might not reach O 
but will settle in L. Similarly if more kinetic energy is 
given to the oscillator, it might reach R instead of O. 

Another important point is that the impulse given to 
the resonator must be very short so that it relaxes fol- 
lowing a free dynamics. The typical frequencies of these 
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FIG. 5: Parameter plane (q, a) showing the difTerent final 
states reached by the system, the left focus * (blue online), 
the center focus + (red online) and the right focus x (green 
online). The parameter /? — 9.76. 



devices are about 500 GHz so the impulse must be around 
5 Thz which is close to the frequency provided by a laser. 
This seems to indicate that an optical pulse generated by 
a laser would be the optimal candidate to switch the de- 
vice. 

To prepare the artificial in a given state, one needs to 
know if this state is really reached. For that one can use 
the pump-probe approach: a first pulse is sent to shift 
the system in the right state, then a small second pulse 
is sent to analyze the state by reflection or transmission. 
This is the object of the next section. 



V. MICROWAVE SPECTROSCOPY OF THE 
SPLIT-RING RESONATOR 

We consider here that the split-ring Josephson res- 
onator is subject to irradiation by a microwave field and 
compute using a scattering theory formalism the response 
of the system. This field could be microwave radiation 
from a wave-guide or it could be a laser beam shining on 
the device. The equations describing the system light- 
ring are the the generalized pendulum equation for the 
fiux ^ and the the Maxwell equation for the electro- 
magnetic field 



V X E 



-Ht-M(, VxH = E( 



(24) 



where E is the electric field, H the magnetic field and M 
the magnetization. Taking the curl of the second equa- 
tion we get the vector wave equation 



V X (V X H) = Htt -f Mf 



(25) 



If we assume that the wave propagates along z and is 
transversely polarized so that H is parallel to x, the nor- 
mal to the plane of the split ring, and E is parallel to y. 
Then we get the scalar wave equation for H 



AH - Hf, 



(26) 



The magnetization is related to the current / in the loop 
by 



M = SI, 



(27) 



where S is the surface enclosed by the ring. Combining 
(PS)) with (P7|) and recalling the expression of the current 
/ given by the term in brackets in ^ we get the final 
system of equations 



LiC<i>tt 



+ Hu I6iz), 



R 



s 



(28) 



where I is the film thickness. 

We introduce the units of flux, magnetic field, time 
and space 



ft 

2^' 



1 



^0 = 



ujt 



(29) 



With these units we normalize time, space, the phase and 
the field as 



H 



H 

is then 



z 1^ $ 
The normalized system obtained from 

HCC - Hrr = -1 (-$rr + Hrr) ^(C) 

$rr + a^r + P sin($) + ^ = H, 
where we have introduced 

lotL „ LIr I 



(30) 



R 



7 



(31) 



(32) 



We assume that the ring is submitted to a fixed mag- 
netic field hg to which it responds with a constant fiux (j)s . 
Then we send in a small electromagnetic pulse 5H and 
examine the response 5$ of the ring using the scattering 
theory. The linearized equations for 5H , (5$ read 



(33) 



Sihrr + Cw5$r /? COs((/)s )(5$ + (5$ = 5H . 
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We now assume periodic solutions 

and obtain the reduced system 

hi^C + ^"^h = (-0 + h) 6{C}, 
[— w'^ + iauj + \ + 13 cos{4>s)\ (j) ~ h. 



(34) 



(35) 



In the scattering we assume the electromagnetic wave to 
be incident from the left of the film located at C = 0. We 
then have 

h = e-'"^^ + i?e*""^ , C < ; Te"*""^, C > , (36) 

where R is the amplitude of the reflected wave and T 
the amplitude of the transmitted wave. We have the 
following interface conditions at C = 



HQ-) - h{Q+), [H]l- = u?-, (-0(0) + hm . 
They imply the two equations for i? and T 

1 + = T, 
-T-{-l + R) = -iujfT 



-u)"^ + iau) + 1 + /3 cos( 



(37) 



+ 1 



from which we obtain the transmission coefficient, 

2(— + 1 + 13 cos 0s + iijja) 



T = 



the reflection coefficient 



D 



R = 



—ct^uP" + i^uj{f3 cos 05 — w^) 



(38) 



(39) 



and where the denominator is 



D = 2(-w^+I+^cos0s)+a7'^^+»[2Q;a;-7a;(/3cos 0s-w^)]. 

(40) 

(41) 



The square of the modulus of R is 

a^7^a;^ + ^"^uj^ {(3 cos 0s 



,2\2 



As seen in section 3, the split-ring oscillator has a finite 
number of equilibria 0s depending on the parameter /3. 
As an example we consider /? = 30 for which the potential 
V{4>) is shown in Fig. [S) The square of the modulus of 
the reflection coefficient (|41l) is plotted in Fig. [7] for the 
five different equilibria. For oj ^ |i?p^Oasa;^, 
for a; — ^ oo — > 1. At some cOg the transmission goes 
to 0, i.e. the medium becomes transparent. Notice the 
difference with a real atom which would absorb incident 
radiation for certain frequencies. The expression for these 
resonant frequencies ujs can be obtained by considering 
the minima of These correspond to the second term 
in the numerator of (|4T|) being zero. We get 



LUs = \/ 13 COS 0s 



(42) 



In the example shown, the spectroscopy data (cjs,0s) is 
given in table 1. 



> 400 
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FIG. 6: Plot of the potential V"(0) for /? = 30. 




FIG. 7: Square of the modulus of the reflection coefficient \R^ 
as a function of the frequency lo for the five different equilibria 
0s = continuous line (red online) and 0s = 27r(l — 1//3) 
dashed line (green onfine). The parameters are /3 = 9.76, 7 = 
1 and Q = 0.01. 



CJs 


5.477 


5.420 


5.238 


4.888 


4.169 


0s 





6.08 


12.15 


18.2 


24.18 



TABLE I: Spectroscopy data (cjs, 0s) for a split-ring oscillator 
with five steady states. The parameters are /3 = 30, 7 = 
1, a = 0.01. 
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VI. CONCLUSION 

We have derived and analyzed a model for split ring 
Josephson resonator or RF SQUID in the superconduct- 
ing regime. If the parameters of the device are chosen 
appropriately, there exist excited states whose number 
can be controlled by carefully tuning the inductance and 
capacity of the ring. We assumed that there are just 
two excited states and showed how an incident magnetic 
flux can shift the system from the ground state to one of 
these excited states. The existence of these excited states 
makes this system similar to an artificial atom with dis- 
crete energy levels. The Josephson oscillator is a unique 
nonlinear element which allows this. Other nonlinear el- 
ements like a diode [l^ , a Kerr material [ll| or a laser 
amplifier would not give these discrete levels. In 
addition, since the oscillator is operating in the super- 
conducting regime the losses are very small as opposed 



to the current meta-materials. 

By sending a microwave field on the resonator we can 
perform a spectroscopy of it and characterize in which 
state it is. Using a scattering theory formalism we com- 
pute the reflection and transmission coefficients for the 
wave. These coefficients differ clearly whether the system 
is in the ground state or in an excited state enabling to 
distinguish them. 
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